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Abstract — The energy storage and attitude control 
subsystems of the typical satellite are presently distinct and 
separate. Energy storage is conventionally provided by 
batteries, either NiCd or NiH, and active attitude control is 
accomplished with control moment gyros (CMGs) or 
reaction wheels. An overall system mass savings can be 
realized if these two subsystems are combined using 
multiple flywheels for simultaneous kinetic energy storage 
and momentum transfer. Several authors have studied the 
control of the flywheels to accomplish this and have 
published simulation results showing the feasibility and 
performance. This paper presents the first experimental 
results showing combined energy storage and momentum 
control about a single axis using two flywheels. 

1. Introduction 

The NASA Glenn Research Center has been developing the 
technology necessary for the use of flywheels on future 
spacecraft for the last several years. Flywheels, as energy 
storage devices, offer several advantages over conventional 
batteries. These include a longer life, greater depth of 
discharge, less complex thermal management and easier 
state of charge measurement (rotor speed). In addition, 
flywheels can also be used as momentum wheels for attitude 
control. This presents the opportunity to combine two 
satellite subsystems, energy storage and attitude control, into 
one, resulting in a lower overall system mass. Several 
authors have studied this possibility and theoretically proved 
the feasibility of the idea [1,2, 3, 4]. These papers are 
generally written from the perspective of the attitude control 
engineers. In this paper, the control is developed from the 
power system and motor control perspective. 

To combine energy storage and attitude control, there must 
at a minimum be one flywheel for each degree of freedom 
required. For example, for one axis of attitude control plus 
energy storage there must be at least two flywheels; for three 
axes of attitude control plus energy storage there must be at 
least four flywheels. More flywheels may be desired in the 
final configuration for redundancy purposes. The control 
laws in any case would be similar — the flywheels must be 


controlled to simultaneously provide the required torque on 
the spacecraft for attitude control and provide the correct 
power to the system as an energy storage device. In addition, 
these functions should be independent, so that a change in 
the power required by the system does not result in a 
disturbance torque to the spacecraft and a change in the 
attitude of the spacecraft does not result in a change in the 
electrical power available to the loads. 

The NASA Glenn Research center has taken the first steps 
toward demonstrating this technology. Successful energy 
storage capability was first experimentally demonstrated 
with a single flywheel in charge, discharge and charge 
reduction modes [5,6,7]. Next, a control algorithm was 
developed and simulated to combine single axis attitude 
control and energy storage [8]. Finally, an experimental 
facility with two flywheels was constructed so that 
simultaneous single axis attitude control and energy storage 
could be demonstrated. Using this facility, this paper will 
present the first experimental verification of the combined 
attitude control and energy storage concept. 

Nomenclature 

Stable Torque applied to air table (N-m) 

T em Electro-mechanical torque produced by motor 

(N-m) 

Fac AC electrical power supplied to motor (watts) 

Vq S Motor q-axis voltage, expressed in the rotor 

reference frame (volts) 

vj s Motor d-axis voltage, expressed in the rotor 

reference frame (volts) 

iq S Motor q-axis current, expressed in the rotor 

reference frame (amps) 

ids Motor d-axis current, expressed in the rotor 

reference frame (amps) 

F qs Motor q-axis inductance (henries) 

F ds Motor d-axis inductance (henries) 
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R s 

CD, 

^af 

^DCtotal 

P 

N 

1> axes 

^flywheels 

Variable* 

ls/a 

hoad 

T* 

1 charge 
Iflywheel 

T D1 

t hss 


Motor resistance (ohms) 

Motor rotor speed (electrical radians/sec.) 

Motor back-EMF constant (volt-sec.) 

Total DC power supplied to both flywheels 
(watts) 

Number of poles of motor 

Desired number of axes of attitude control 

Number of flywheels in a combined attitude 
control/energy storage system 

Commanded value of variable 

Current from solar array or power supply (amps) 

Current into DC load (amps) 

Commanded flywheel charging current (amps) 

Total DC current into flywheel units (amps) 

DC Current into D1 flywheel unit (amps) 

DC Current into HSS flywheel unit (amps) 



Figure l.-Test Facility Schematic 

Electrically, the motor/generator of each flywheel is 
connected to the same dc bus through its own inverter and 
filter as shown in Figure 2. In charge mode, the DC current, 
Iflywheel? is positive and the speed of the flywheels is 
increasing. In discharge mode, the flywheels are decreasing 
in speed and providing power to the DC bus. 


2. Control Algorithm Development 

The inner-most loop of the overall control algorithm is 
based on a motor control technique known as field 
orientation control. This technique allows accurate, high 
bandwidth torque control of the flywheel motor/generator by 
regulating rotor reference frame currents that are DC in the 
steady state. A general description of field orientation of 
permanent magnet machines can be found in [9] and a 
specific discussion of the motor control techniques used in 
the flywheel program can be found in [10,11]. The 
description of the overall flywheel system control given in 
this paper assumes a basic familiarity with the field 
orientation technique. 

Overall Control Derivation 

The single axis system to be controlled is shown 
schematically in Figure 1. It consists of two flywheels 
mounted on an air-bearing table. Table movement is 
restrained by a load cell that is connected between the outer 
edge of the table and a fixed point. The load cell measures 
the net force applied to the table due to the two flywheels. 
The axial torque on the table can then be computed from 
(r x F). The net steady state torque applied to the table is the 
sum of the flywheel motor reaction torques as shown in 
equation 1. 



Figure 2. — Electrical Schematic 

The AC electrical power supplied to the motor can be 
expressed using the motor currents and voltages expressed 
in the rotor reference frame. 

3 

P AC — 2 (vqsiqs Vdslds) (2) 

The motor currents are controlled such that i r ds is regulated to 
zero and thus (2) reduces to (3). 

3 

P AC — 2 (Vqsiqs) (3) 


Ttable “Teml T em 2 


( 1 ) 


Equation (4) gives the motor voltage, Vq S , from [9]. 
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Assuming steady state conditions (under field orientation 
control, the quantity iq S is DC in steady state) and neglecting 
the resistive drop (which is generally small compared to the 
back EMF term, co,A a f), (4) reduces to (5). 

Vqs ~ ®r^af (5) 

Substituting (5) into (3) results in the expression for the 
steady state AC electrical power of the motor/generator. 

* 3. r „ 

P AC ~ 2 iqs 0\^af (6) 

Note that if iq S is positive and the machine is spinning in the 
positive direction as defined by Flywheel 1 in Figure 1, then 
the power is positive. This means that the machine is 
motoring (charging). Similarly, if iq S is negative, the power is 
negative, implying that the machine is generating 
(discharging). Thus charge and discharge of the flywheel 
can be regulated by properly controlling iq S . 

If the inverter/filter losses are neglected, the DC power 
supplied to each flywheel is approximately equal to the AC 
power shown in (6) in steady state. The total flywheel power 
for the system is the sum of the power for each individual 
flywheel. 


P DCtotal ~ P AC1 + P AC2 “ 

3 . r 3 . r 

2 iqsl Orl^afl + 2 ^qs2 ^r2^af2 (7) 

Equations (1) and (7) form the fundamental basis for the 
control of the two flywheels to achieve simultaneous attitude 
control and energy storage. In (1), the control variable is the 
motor/generator torque, x cm , while in (7) the control variable 
is the motor/generator current, iq S . To combine the two 
equations, the relationship between the motor/generator 
torque and current must be used and that is given in (8). 

3 p 

^em _ 2 2 C^ds “ Lq S ) id s ] (8) 

With i^ regulated to zero, it becomes 



-1 -1 



Teml 
_ T e m2_ 




" 3PAafl 

3P 2 X af2 l 


Ttable 


' 4 

' 4 

"iqsl" 

_ P&c- 


2 ®Aafl 

2 ®Aaf2 

_lqs2_ 


( 11 ) 


Note that (10) can be modified to include additional 
flywheels and axes of control by determining the 
relationship between the torque produced by each flywheel 
motor and its effect on each axis. This will depend on the 
geometry of a particular configuration but in general, the left 
hand side vector of (10) will be of length (N axes +1) and the 
right hand side vector will be of length (N flywhee i s ). Equation 
(11) can then be determined from (10) and (9) for an 
arbitrary number of flywheels. 


Determining Commanded Control Quantities 
In our system, fast motor torque control is achieved through 
the use of the field orientation technique and a high 
bandwidth synchronous frame current regulator [12]. Thus it 
is most useful to solve (1 1) for the necessary motor currents 
given a commanded value of table torque and DC power. 
This is shown in (12). 



2 

2CG2^af2 P 2 Af2 

- * “ 

Ttable 

3^ a n^af2(C0iP2-0)2Pl) 

_ -2(Oi?l a fi -P A a fl_ 

LPdcJ 


Equation (12) is the same result first given in reference [8]; 
however, here it is expressed more generally in matrix form. 

The commanded values of x ta bi e and P dc can be the result of 
other control loops or can be commanded directly. For 
example, the table torque command could be the result of a 
position feedback loop to control the position of the table. 
Simulation results for this case are given in [8]. In the results 
presented herein, x ta bie is commanded directly. There is no 
closed loop control on the torque applied to the table. 

The power command, P dc , is the result of the DC bus 
regulation algorithm that is described in detail in [5,6,7]. It 
will also be summarized here. 


Tem — 2 2 (9) 

Combining (1), (7) and (9), the system control can either be 
expressed in terms of motor/generator torque, as in (10), or 
motor/generator current, as in (1 1). 


Power Regulation 

There are three modes of operation for the batteries in a 
spacecraft power system: charge, charge reduction and 
discharge. The control for the flywheel system in this 
experimental demonstration is constructed in the same 
fashion so as to facilitate the eventual demonstration of 
flywheel energy storage on a spacecraft. 
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In charge mode, the solar array produces enough current to 
charge the flywheels at a setpoint and to provide the 
required load current. The solar array electronics regulate 
the DC bus voltage during charge mode. In charge reduction 
mode, the solar array continues to provide load current but it 
cannot provide enough current to charge the flywheels at the 
setpoint. When this occurs, the DC bus voltage regulation 
function is transferred to the flywheel system. Finally, in 
discharge mode, the flywheel system provides all of the load 
current and regulates the DC bus voltage. The features of 
each mode are given in Table 1 using the variables defined 
in Fig. 2 where I charge is the commanded charging current. 


Table 1. — Flywheel System Operating Mode Characteristics 


Mode 

Current 

DC Bus 
Voltage 

Full Sun 
“Charge” 

Is/a = hoad + Icharge 
Ifiywheei — Icharge 

Regulated by 
solar array 
system 

Partial Sun 
“Charge 
Reduction” 

hoad + Icharge > h/a > hoad 
Icharge > Ifiywheei > 0 

Regulated by 
flywheel system 

Eclipse 

“Discharge” 

hoad = " Ifiywheei 
Ifiywheei < 0 

Regulated by 
flywheel system 


In charge mode, the power command, P dc , is determined 
from the combination of a feed-forward term and a feedback 
term as shown in Figure 3. The feed-forward term is simply 
the commanded charging current. In an ideal world, where 
all parameters are known perfectly, this would be sufficient. 
The feedback term is derived from the comparison of the 
measured current and the commanded current. It provides 
error compensation that forces Ifiywheei to track the 
commanded value. Using the feed-forward term allows the 
gains on the PI portion to be set much lower than otherwise 
necessary and results in fast response with less noise on P dc . 



Figure 3. — Charge Mode Controller 

In charge reduction and discharge modes, the flywheel 
system must regulate the DC bus voltage. The block diagram 
for this algorithm is shown in Fig. 4. Again, there are two 
components to the controller: the PI portion and a 
disturbance decoupling (DD) portion. In the decoupling 
portion, the DC current is measured and used as an early 
indicator to the controller whether there has been an increase 
or decrease in load. If there is a sudden increase in load, the 


capacitor will initially maintain the bus voltage and there 
will be an increase in the DC current, Ifiywheei? to supply the 
new load. This increase in Ifiywheei is used by the controller 
increase the DC power command, P dc , even before a drop in 
the DC bus voltage causes the PI portion to respond. 



Figure 4. — Discharge and Charge Reduction Controller 


In the PI portion the measured DC bus voltage is compared 
to the commanded DC bus voltage and the error signal is 
inverted. This is because (using the direction references of 
Fig. 2) if the DC bus voltage is decreasing, there needs to be 
an increase in the average inverter current to charge the 
capacitor and increase the bus voltage, and if the DC bus 
voltage is increasing, there needs to be a decrease in the 
inverter current. The PI portion makes up for any 
inaccuracies in the disturbance decoupling portion and 
maintains the DC bus voltage at the set point. 

The transition from current regulation (Fig. 3) to voltage 
regulation (Fig. 4) is accomplished in the following manner 
and shown in Fig. 5. The solar array regulates the bus 
voltage to a set point value higher than the flywheel 
regulation set point as long as the solar array current is 
sufficient to provide both the load and the charging current, 
I charge • Once the solar array current begins to drop off, the 
DC bus voltage begins to fall and the flywheel current, 
Ifiywheei? also drops. This transition is detected in the 
controller by comparing the difference between the actual 
DC bus voltage and the flywheel set point voltage to the 
’’voltage transition constant”, VTC, as seen in Fig. 5. Once 
this difference is less than the VTC, the integrator in the PI 
portion of the controller is reset. This reduces the if nv 
command at point 2 to a value slightly larger than Ifiywheei- 
This value is then compared to the charge current set point, 
I charge • If it is l es s than I charge ? which it will be if the solar 
array is not producing enough current, then the system 
transitions into charge reduction mode where the DC bus 
voltage is regulated by the flywheel system. 

Similarly, as the system moves from eclipse into sunlight, 
the solar array will produce more and more current. When 
the solar array produces enough current to meet the load 
demand, the i f nv command at point 2 in the controller will 
become positive. When it exceeds the charge current set 
point, I charge ? the integrator in the current regulator portion 
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of the controller is reset and the system transitions back into 
charge mode where the flywheel system regulates the current 
into the flywheel and the solar array system regulates the DC 
bus voltage. 



Figure 5. — Combined charge/discharge power 
regulation control block diagram. 


It is important to note that the three modes of operation, 
charge, charge reduction and discharge, were originally 
defined based on a battery energy storage system. The 
flywheel energy storage system is capable of regulating the 
DC bus voltage at all times, both when charging and 
discharging, obviating the need for current and voltage 
regulation modes and the transition between them. This 
would result in an overall simpler control strategy, even 
when considering the necessary provisions to prevent over- 
speed or over-current operation. This modification of the 
control will be investigated in future efforts. 


Overall Control Block Diagram 

The overall system control block diagram is given in Figure 
6. Each portion of the controller has been clearly defined 
and labeled, with references to the corresponding equations 
if appropriate. Although complex, the control can be broken 
down into manageable blocks as shown such that 
improvements in a particular area can be individually 
developed. For example, a different power and attitude 
control algorithm block could be evaluated while the rest of 
the control remained the same. 

3. Simulation and Experimental 
Assumptions and Conditions 

The next four sections will present a combination of simulation 
and experimental results. First, the simulation and experimental 
conditions and equipment will be described in this section. In 
the next section, the baseline set of simulation and experimental 
results will be given for the control as described thus far, in 
particular, Fig. 6 and equation (12). It will be seen that there are 
some transient disturbances that occur on the power bus due to a 
step change in the table torque command. Next, in Section 5, 
simulation and experimental results will be presented for the 
case where the control algorithm is modified to include the 
motor and AC side resistances. This modification improves the 
disturbance rejection of the power system bus during table 
torque command transients. Finally, in Section 6, a set of 
experimental results will be presented to show the effect of 
spinning the two flywheels “backwards” relative to the first two 
sets of data. 
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The simulations were done using Matlab/Simulink and were 
based on the block diagram of Figure 6. The PM motors 
were simulated based on standard motor equations [8]. The 
switching of the voltage source inverter was not modeled 
(PWM harmonics were not considered) and the stator 
resistance of the motor model was increased to account for 
the switch, conductor and filter resistances for the simulation 
results in Section 5. 

The important control parameters for the two flywheels, 
known as High Speed Shaft (HSS) and Development 1 (Dl), 
are given in Table 2. It can be seen that the two flywheels 
are very dissimilar in terms of voltage rating, inertia, back 
EMF constant and number of motor poles. The use of these 
two flywheel units was dictated by budget and availability 
constraints. In any actual flight system, the two flywheels 
would be designed identically and any differences would be 
due only to manufacturing tolerances. Although less than 
ideal, the fact that we can control the system with two 
different motors and two very different inertias is a 
testament to the robustness of the control technique. 


Table 2. — Parameters of Flywheel Motors 


Dl Motor/Generator Parameters 

J 

Rotor polar moment of inertia 

0.0664 kg-m 2 

L q 

Motor q-axis inductance 

25 pH 

L d 

Motor d-axis inductance 

19 pH 

R* 

Stator resistance 

0.02 Q 

^af 

Back EMF constant 

0.0103 volt-sec. 

P 

No. of poles 

2 

V n om 

Rated L-N RMS voltage 

45.7 volts 

High Speed Shaft (HSS) Motor/Generator 
Parameters 

J 

Rotor polar moment of inertia 

0.00377 kg-m 2 

L q 

Motor q-axis inductance 

142 pH 

L d 

Motor d-axis inductance 

101 pH 

R* 

Stator resistance 

0.035 Q 

^af 

Back EMF constant 

0.0144 volt-sec. 

P 

No. of poles 

4 

V n om 

Rated L-N RMS voltage 

128 volts 


The experimental tests were limited to speeds of 12,000 rpm 
or less because the magnetic bearings for the system had not 
yet been optimally tuned for higher speed operation. The 
initial simulation speeds were set to 1 1,000 rpm to be in the 
range of the experimental conditions. 



The air table dynamics, with the load cell attached, were 
measured and then modeled to fit the measured data. The air 
table was moved from an equilibrium position and allowed 
to return to equilibrium without the influence of the flywheel 
motors as seen in Fig. 7. This response was modeled as a 
second order system with a damping factor of .0186 and a 
natural frequency of 8.57 Hz. Using these factors, the 
simulated response is also shown on Fig. 7. 

Each set of data is grouped according to the three general 
test conditions listed in Table 3. These tests demonstrate a 
representative range of functionality for the combined 
flywheel system energy storage and single axis attitude 
control. For each test, a set of data is presented that includes 
the DC currents (defined in Fig. 2), the DC bus voltage, the 
torque on the air table, the motor speeds and the motor 
phase currents. 


Table 3. — General Test Conditions for Simulation 
and Experimental Results 


Test 

Number 

Power 

Regulation 

Mode 

Table Torque 
Command (N-m) 

Load 

1 

Charge to 
Discharge 

T* = 0 

Step 

change 

2 

Charge 

T* = 0 ^ +0.5 — » 0 

Constant 

load 

3 

Discharge 

T* = 0 — > -0.5 — > 0 

Constant 

load 


4. Simulation and Experimental Results: 
Part A 


Test 1 

In this section, the functionality of the control in the charge 
and discharge modes with a zero torque command on the 
table is shown with simulation and experimental results. 
Figures 8 through 13 give the simulation results and Figures 
14 through 19 give the experimental results for the Test 1 
conditions. Point “a” is a transition from one commanded 
current value to another while in charge mode, point “b” is 
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a transition from charge to discharge mode, point “c” is a 
step change in the electrical load and point “d” is a transition 
from discharge to charge mode. (Note: point “d” is not 
shown in the experimental results because the time capture 
setting was too small. This will be modified in the next 
publication on this work.) The table torque command is set 
to 0 N-m throughout these tests. 

Figure 8 shows the total flywheel current following the step 
change in the charging current command. It also shows the 
source current transitioning to zero and the flywheel current 
going negative to supply the entire load as the system moves 
from charge to discharge mode. A negative system flywheel 
current indicates that the flywheels are providing the power 
to the system; a positive flywheel current indicates that 
energy is being stored in the flywheels. Figure 14 gives the 
same information for the experimental case. 

Figures 9 and 15 show the split in the total flywheel system 
current between the two flywheels, HSS and Dl. In the 
experimental case (Figure 15) the HSS and Dl currents are 
practically equal so that one masks the other in the plot. 
Figure 10 shows the voltage change from the charge value of 
125 volts to the flywheel regulated discharge value of 120 
volts. It can be seen that there is a slight voltage transient 
when the power load is increased at point “c”. Figure 15 is 
the corresponding figure for the experimental case. 

Figure 1 1 shows the flywheel rotor speeds. It can be seen 
that the HSS flywheel decelerates much more quickly than 
the Dl flywheel to provide the required power to the load 
while maintaining 0 torque on the table as seen in Figure 12. 
Figures 17 and 18 are the corresponding experimental 
results. 

Finally, Figure 13 shows the Dl motor phase current 
envelope for simulation and Figure 19 is the corresponding 
experimental plot. These plots demonstrate the functionality 
of the control in the charge and discharge modes with a zero 
torque command on the table. 

Test 2 

The next set of data demonstrates the functionality of the 
control in charge mode with a step change in the table 
torque. The Test 2 simulation results are presented in 
Figures 20 through 25 and experimental results are given in 
Figures 26 through 31. A step change in the charging current 
command is given at point “a”, then a torque step command 
of +0.5 N-m is given at point “b”, and at point “c” the torque 
command is returned to 0. (The step change in current 
command at point “a” is not shown in the experimental 
results.) 

Figure 20 shows the total flywheel current in charge mode 
during a step change in the table torque command. Figure 26 
presents the corresponding experimental results. In both 
plots, it can be seen that a step change in the table torque 
results in a disturbance to the flywheel charging current. 


The action of the PI control in the charge mode portion of 
the power controller (Fig. 3) is seen to bring the flywheel 
current back toward the commanded value. This disturbance 
is caused by an inaccurate P dc command to the attitude 
control and bus regulation algorithm. The control algorithm 
is modified to correct this in the results presented in the 
Section 5. 

Figures 21 and 27 show the individual flywheel currents 
during the torque step command for simulation and 
experimental results respectively. It can be seen that the 
HSS is discharging (negative DC current) while the Dl unit 
is charging (positive DC current) to exchange the correct 
amount of momentum to apply the commanded torque to the 
table. Figures 23 and 29 also show the exchange of 
momentum between the HSS and Dl because HSS is 
decelerating (discharging) while Dl is accelerating 
(charging). 

Figures 24 and 30 show the table torque response to the step 
change command. Both simulation and experimental results 
show the torque following the commanded value with the 
expected dynamic response of the air table and load cell. 

Finally, figures 25 and 31 show the Dl motor phase current 
during the transitions. 

Test 3 

This set of data demonstrates the functionality of the control 
in discharge mode with a step change in torque. To 
demonstrate torque control in both the positive and negative 
directions, a negative torque command is given for this case. 
The Test 3 results are presented in Figures 32 through 37 for 
simulation and Figures 38 through 43 for the experimental 
results. Point “a” corresponds to a step change in the 
charging current command, point “b” is the transition from 
charge mode to discharge mode, and point “c” is the torque 
step from 0 to -0.5 N-m. (Again, the step change in current 
command at point “a” is not shown in the experimental 
results.) 

Figures 32 and 38 present the DC currents for simulation 
and experimental results respectively. The total flywheel 
current becomes negative when the system moves into 
discharge mode. Figures 33 and 39 show the individual 
flywheel currents during the torque step command for 
simulation and experimental results respectively. The HSS 
DC current is positive for a negative torque step and the Dl 
current is negative indicating that the HSS flywheel is 
charging while the Dl flywheel is discharging. 

Figures 35 and 41 show the corresponding flywheel speeds 
for simulation and experimental results respectively. Both 
flywheels initially start out discharging, but when the torque 
step is commanded, the HSS accelerates while the Dl unit 
discharges more quickly to provide both the electrical power 
and momentum transfer to the HSS necessary to meet the - 
0.5 Nm torque command. 
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Figure 8. — Source, Load and Total Flywheel 
System DC Currents for Test 1 Simulation 



Figure 10. — DC Bus Voltage for Test 1 Simulation 
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Figure 12. — Table Torque for Test 1 Simulation 



Figure 9. — D1 Flywheel, HSS Flywheel and Total Flywheel 
System DC Currents for Test 1 Simulation 



Figure 11. — Flywheel Rotor Speeds for Test 1 Simulation 
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Figure 13. — D1 Motor Phase Current for Test 1 Simulation 
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Figure 14. — Source, Load and Total Flywheel System 
DC Currents for Test 1 Experimental Results 
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Figure 16. — DC Bus Voltage for Test 1 Experimental Results 
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Figure 18. — Table Torque for Test 1 Experimental Results 
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Figure 15. — D1 Flywheel, HSS Flywheel and Total Flywheel 
System DC Currents for Test 1 Experimental Results 



Figure 17. — Flywheel Rotor Speeds for Test 1 
Experimental Results 
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Figure 19. — D1 Motor Phase Current for Test 1 
Experimental Results 
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Figure 20. — Source, Load and Total Flywheel Figure 21. — D1 Flywheel, HSS Flywheel and Total Flywheel 

System DC Currents for Test 2 Simulation System DC Currents for Test 2 Simulation 



Figure 22. — DC Bus Voltage for Test 2 Simulation 



Figure 23. — Flywheel Rotor Speeds for Test 2 Simulation 



Figure 24. — Table Torque for Test 2 Simulation 
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Figure 25. — D1 Motor Phase Current for Test 2 Simulation 
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Figure 26. — Source, Load and Total Flywheel System DC 
Currents for Test 2 Experimental Results 


Figure 27. — D1 Flywheel, HSS Flywheel and Total Flywheel 
System DC Currents for Test 2 Experimental Results 




Figure 28. — DC Bus Voltage for Test 2 
Experimental Results 


Figure 29. — Flywheel Rotor Speeds for Test 2 
Experimental Results 




Figure 30. — Table Torque for Test 2 Experimental Results Figure 31. — D1 Motor Phase Current for Test 2 

Experimental Results 
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torque, N-m UU bus voltage, volts 



Figure 32. — Source, Load and Total Flywheel System DC 
Currents for Test 3 Simulation 



Figure 33. — D1 Flywheel, HSS Flywheel and Total Flywheel 
System DC Currents for Test 3 Simulation 



Figure 34. — DC Bus Voltage for Test 3 Simulation 



Figure 35. — Flywheel Rotor Speeds for Test 3 Simulation 
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Figure 36. — Table Torque for Test 3 Simulation 
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Figure 37. — D1 Motor Phase Current for Test 3 Simulation 
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Figure 38. — Source, Load and Total Flywheel System DC 
Currents for Test 3 Experimental Results 


Figure 39. — D1 Flywheel, HSS Flywheel and Total Flywheel 
System DC Currents for Test 3 Experimental Results 



Figure 40. — DC Bus Voltage for Test 2 Experimental Results 


Figure 41. — Flywheel Rotor Speeds for Test 2 
Experimental Results 
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Figure 42. — Table Torque for Test 3 Experimental Results Figure 43. — D1 Motor Phase Current for Test 3 

Experimental Results 
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Figures 34 and 40 show the DC bus voltage. The torque step 
occurs while the system is in discharge mode (the flywheel 
system is regulating the DC bus) and it can be seen that there 
is a disturbance on the bus voltage at the torque step. Similar 
to the Test 2 results, this occurs because an inaccurate P d * 
command is given to the attitude control and bus regulation 
algorithm. This also is corrected by the control algorithm 
change presented in the next section. 


(discharge mode) during a step change in the motor torque. 
The revised equation for the AC power to include the 
resistive drop is given in (13). 

3 

P AC ~ 2 (Or^af - * - iqs Rs) (13) 

The approximate total DC power then becomes 


Figures 36 and 42 show the table torque response to the step 
change command. In the simulation case, the torque follows 
the commanded value exactly. However, experimentally, a 
steady state error is observed. This could be caused by using 
the wrong motor parameters in the control equations, 
specifically the back EMF constants A, afm and ?t afHSS could 
be somewhat inaccurate. This is investigated in Section 6. 

Figures 37 and 43 show the D1 motor phase current 
envelope during the transitions. 

This initial series of tests, demonstrated both in simulation 
and experimentally, shows the functionality of the controller 
through several different combinations of power modes and 
torque steps. Three discrepancies were noted: there was a 
disturbance in the charging current during charge mode 
when a torque step was applied, there was a disturbance in 
the bus voltage during discharge mode when a torque step 
was applied and the table torque did not reach the 
commanded negative value in the experimental results. The 
first two of these discrepancies, power disturbances due to 
torque steps, are addressed in the next section. 

5. Simulation and Experimental Results: 
PartB 

In this section, the effect of modifying the control algorithm 
to include the AC resistance neglected in (4) is investigated. 
It will be seen in both simulation and experimental results 
that the modification of the control algorithm improves the 
power bus disturbance rejection during torque steps in both 
charge and discharge modes of operation. 

Control Modification 

In the original derivation of the control, the steady state 
motor voltage was calculated from (5) which neglected the 
AC voltage drop across the motor resistance. This was 
because it is generally true that this drop is small relative to 
the back EMF term, Civ^af- 

Vqs ~ COr^af ( 3 ) 

However, in this case, the iq S R s term represents not only the 
voltage drop across the motor/generator resistance, but also 
the drop across the inverter switches, the AC filter and 
conductors leading to the machine. Neglecting it in this case 
turned out to unexpectedly lead to a disturbance in either the 
charging current (charge mode) or DC bus voltage 


DC total ~ PaCI + P AC2 “ 

3 3 

2 Iqsl^CVl^afl^iqslRsl) 2 iqs2(c0r2^af2 _ ^iqs2Rs2)- (14) 


This expression is then combined with (1) and (7) to form 
the relationship between the DC power, the table torque and 
the motor currents as shown in (15). 


3 PiA, a fi 

4 



3P2^af2 

4 


3 3 

2((0 r i^afl _ ^iqsl Rsl) 2C^r2^af2~^iqs2 Rs 2 ) 


I” Iqsl 
_lqs 2 _ 


(15) 


It can be seen that the inclusion of the resistance term leads 
to a nonlinear expression relating the control variables, the 
motor currents, to the desired values of table torque and DC 
bus power. Thus it becomes more difficult to determine the 
commanded motor current values from the desired values of 
table torque and DC bus power. One approach to solving 
this is to linearize (14) about an operating point defined by 
the measured values of currents and then use the linearized 
equation in conjunction with equations (1) and (7) to form 
the matrix expression. 


However, the approach used in this case is to use the matrix 
expression given in (15) with the modification shown below. 
The iq S values used in the matrix now represent the measured 
(i q r s M ) values fed back from the current sensors on the motor 
phases and the iq S values used in the vector are the 
commanded (iq S ) ones. This is shown in (16). 

Ttable ___ 

_ Pdc_ “ 

3Pl^ a fl 3P2?laf2 

4 " 4 

3 3 

2(c0 r i^ a fi+iqsi Rsl) 2C^r2^af2~^iqs2 Rs2) 

(16) 

Equation (16) is solved for the commanded currents as 
shown in (17). The result given in (17) is almost identical to 
the linearized solution, however, it is easier and more 
intuitive to implement. 
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Change in Table Torque with Original Control 
(eq. (12)), Simulation Results 



Figure 45. — Charging Current for Step Change in Table 
Torque with Modified Control (eq. (17)), 
Simulation Results 



Figure 46. — Disturbance in Charging Current for Step 
Change in Table Torque with Original Control 
(eq. (12)), Experimental Results 


Figure 47. — Charging Current for Step Change in Table 
Torque with Modified Control (eq. (17)), 
Experimental Results 
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(eq. (12)), Simulation Results 
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Figure 50. — Disturbance in DC Bus Voltage for Step Change 
in Table Torque with Original Control (eq. (12)), 
Experimental Results 



Figure 49. — DC Bus Voltage for Step Change in Table 
Torque with Modified Control (eq. (17)), 
Simulation Results 
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Figure 51. — DC Bus Voltage for Step Change in Table 
Torque with Modified Control (eq. (17)), 
Experimental Results 
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3[P2^af2((Ol^afl+ IqslRsl) " Pl^afl(t02^af2+iqs2Rs2)] 


X 


2(C02^af2+iqs2Rs2) P 2 Vil2 

-2 (COl A-af 1 +lqsl Rs 1 ) -PlA-afl. 


r * 

Stable 

LPdcJ 


(17) 


direction and the HSS flywheel was spinning in the counter- 
clockwise (CCW) direction. The torque step experiment as 
repeated for D1 spinning in the CCW direction and HSS spinning 
in the CW direction. If the steady state error was due to a control 
algorithm problem or inaccurate estimate of the motor parameters, 
changing the spin direction should result in a steady state error for 
a positive torque step and no error for a negative step. That was 
not the case, as seen in Figures 52 through 59. 


To verify the improvement in the control by using the modified 
equations, the two torque step tests (Test 2 and 3) were repeated 
with the adjusted control. Test 2 is the case where the torque step 
occurs while in charge mode. Simulation and experimental results 
for the charging current during the torque step are shown in Figs. 
44 through 47 with the previous results (without the control 
modification) shown also for comparison. It can be seen that the 
modification reduces the charging current disturbance significantly 
both in simulation and experimentally. 

Test 3 is the case where the torque step occurs while the system is 
in discharge mode. Simulation and experimental results for the DC 
bus voltage during the torque step are shown in Figs. 48 through 
51, again with the original results repeated for comparison. It can 
be seen that the modification reduces the DC bus voltage 
disturbance significantly during the torque step transient. 

6. Experimental Results: Part C 

In this section, the cause of the steady state torque error observed 
experimentally for a negative step change in torque (Test 3, Fig. 
42) is investigated. It is concluded that the error is not due to an 
error in the control algorithm or estimated motor parameters but 
rather is due to the mechanical characteristics of the load cell being 
used to measure the table torque. 

To investigate the steady state torque error, Test 3 was conducted 
again (experimentally only because the error was not observed in 
simulation) but with each flywheel spinning in the opposite 
direction from the previous tests. In the Part A and Part B results 
presented in Sections 4 and 5, the D1 flywheel was tested while 
spinning in the clockwise (CW) 


Figures 52 through 55 give the results for a positive step change in 
torque with the HSS spinning CW for Fig. 52 and 54 and CCW 
for Fig. 53 and 55. It can be seen that the positive step in torque 
exhibits no steady state error regardless of the direction of spin and 
whether the flywheels are accelerating or decelerating during the 
step. 

Figures 56 through 59 give the results for a negative step change 
in torque with the HSS spinning CW for Fig. 56 and 58 and CCW 
for Fig. 57 and 59. In this case, there is a steady state error in the 
torque for either spin direction. The conclusion is that this error is 
due to the mechanical characteristics of the load cell and not due 
to any control inaccuracy. 

7. Conclusions 

The NASA Glenn Research Center has experimentally 
demonstrated the feasibility of combined single axis attitude 
control and energy storage using two flywheels. Although the 
concept of combining attitude control and energy storage has been 
studied previously in the literature, this is the first experimental 
demonstration of such a system. Torque steps during charge and 
discharge modes and load (power) steps while maintaining DC 
bus regulation have all been demonstrated. Additionally, 
decoupled torque and power control was improved by modifying 
the original control to include the effects of resistive drops on the 
AC side of the system. The overall system control is accomplished 
by the integration of three techniques: a position sensorless field 
orientation algorithm for the inner loop motor control, a 
charge/discharge algorithm for power regulation and an outer loop 
control to simultaneously achieve commanded single axis torque 
and electrical power. 
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Figure 52. — Positive Torque Step for CW Rotation of HSS 
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Figure 53. — Positive Torque Step for CCW Rotation of HSS 
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Figure 54. — HSS Accelerating During Torque Step Figure 55. — HSS Decelerating During Torque Step 
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Figure 56. — Negative Torque Step for CW Rotation of HSS Figure 57. — Negative Torque Step for CCW Rotation of HSS 



Figure 58. — HSS Decelerating During Torque Step 



Figure 59. — HSS Accelerating During Torque Step 
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